Long-range plasmon-assisted energy transfer over doped graphene 



Kirill A. VelizhanirH 

Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA 

V. Shahbazyar|t] 

Department of Physics, Jackson State University, Jackson, Mississippi 39217, USA 

We demonstrate that longitudinal plasmons in doped monolayer graphene can mediate highly 
efficient long-range energy transfer between nearby fluorophores, e.g., semiconductor quantum dots. 
We derive a simple analytical expression for the energy transfer efficiency that incorporates all the 
essential processes involved. We perform numerical calculations of the transfer efficiency for a pair 
of PbSe quantum dots near graphene for inter-ffuorophore distances of up to 1 /xm and find that 
the plasmon-assisted long-range energy transfer can be enhanced by up to a factor of ^lO'^ relative 
to the Forster's transfer in vacuum. 



I. INTRODUCTION 

Forster resonance energy transfer (FRET)^ between 
spatially separated donor and acceptor fluorophores, 
such as dye molecules or semiconductors quantum dots 
(QD), underpins diverse phenomena in physics, chem- 
istry and biology. Examples include photosynthesis, ex- 
citon transfer in molecular aggregates, interactions be- 
tween proteins^ ^ and, more recently, energy transfer be- 
tween QDs and QD-protein assemblies .^^^ During the 
past decade, remarkable progress has been made in appli- 
cations of FRET spectrosco py, e .g., in protein folding,''!^ 
live cell protein localization,'^^ biosensin^^^^^ and light 
harvesting. The range of present and potential appli- 
cations of FRET is, however, limited by its intrinsi- 
cally short-range nature. Indeed, the underlying FRET 
mechanism - the direct Coulomb interaction between flu- 
orophores - supports efficient transfer only at donor- 
acceptor distances {vad) below the typical Forster ra- 
dius oi rp ~ 10 nm.^ At larger distances, the Coulomb 
potential between electrically neutral donor and accep- 
tor decreases rapidly, and the FRET efficiency falls off 
as ^T^/r^d' Substantial efforts have been undertaken 
to improve the efficiency and increase the range of en- 
ergy transfer (ET) at the nanoscale by utilizing surface 
plasmons (SP) and surface plasmon-polaritons (SPP) 
as inter mediaries.'^^t^ Placing molecules or QDs near 
a metal film or a nanoparticle can lead to a signifi- 
cant improvement of ET efficiency (ETE) - the frac- 
tion of donor's energy transferred to the acceptor .'^I^'^ 
In metals, however, the efficiency of plasmon-mediated 
ET channels is limited by significant Ohmic losses and 
plasmon-enhanced radiative losse^^ result ing in a rela- 
tively modest (-10) overall ETE increas^O™^ or 
even its reductiortf^SlI near metal structures. 

In this article, we propose to exploit collective ex- 
citations in graphene as efficient ET intermediaries. 
Graphene has recently emerged as a novel intrinsically 
two-dimensional material^^ with unique electronic and 
optical properties. Clean graphene samples are char- 
acterized by long electron scattering times and much 
lower, compared to metals, Ohmic losses due to relatively 
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FIG. 1. (a) Schematic of ET between donor (D) and acceptor 
(A) through the GP excitation, (b) Density plot of the imag- 
inary part of graphene density correlation function (darker 
shades of gray correspond to higher magnitudes). The Fermi 
level is set to e/ = 0.6 eV. (c) Normalized donor energy loss 
rate versus donor-graphene distance for suspended graphene 
in vacuum. The energy loss rate is plotted for three different 
excitation energies. Thick solid portions of lines mark the in- 
tervals of the exponential decay of kd - intervals where GP 
is most efficiently excited (tagged "plasmon" in the graph). 
The inset shows the zoomed-in left-top portion of the main 
graph in the panel (c). Intervals of efficient electron- hole pair 
excitation at low Zd are tagged "e-h pairs" in the inset. 



weak electron-phonon interaction.^^ Furthermore, doped 
graphene possesses a stable in-plane longitudinal plas- 
mon in the infrared frequency range with gate-tunable 
wavelength, Ap, well below radiation (or SPP) wave- 
length Ao at the same frequency.'^ Recent optical imaging 
of graphene plasmons (GP) propagating in a graphene 



2 



ribbon on SiC substrate indeed demonstrated the high 
degree of GP locahzation char acteri zed by hght-to-GP 
wavelength ratio of Aq/A^ ^ 40.'^^l^ The large GP local 
density of states, as compared to that of SPP, permits 
very efficient GP excitation by a local probe such as a 
scanning tunneling microscope (or atomic force micro- 
scope) tip, an excited molecule or QD placed at a close 
distance {^Xp) from the graphene sheet. ^^ ^^ Recently, 
superradiance from two emitters placed near graphene 
was studied, and it was demonstrated that the inter- 
action of fluorophores with plasmons in graphene can be 
strong enough to significantly enhance the superradiant 
coupling between these fluorophores. 

In this work, we demonstrate that plasmons in doped 
graphene can mediate a highly efficient long-range ET 
between fluorophores, e.g., QDs. A photoexcited donor 
with energy hu situated at a distance Zd from the 
graphene sheet excites a GP which propagates a distance 
Rad ^ Ap in the plane before exciting a remote acceptor 
at a distance Za from graphene (see schematics in Figurejl] 
and Figure [2|. Importantly, when the GP wave reaches 
the acceptor, its intensity is reduced only by factor of 
(xXp/Rad due to the strictly in-plane GP propagation. 
This, along with the efficient fluorophore-GP coupling at 
^o,d ^ Ap, leads to a very strong ET enhancement (up to 
~10^) as compared to the FRET channel, at distances 
far exceeding the Forster radius. 

We show that at large transfer distances, Rad ^ Ap, 
ETE between donor and acceptor is given by (see Sec- 
tion |ll] for the detailed derivation) 
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is the characteristic ET length which, in high mobility 
graphene, only weakly depends on Rad- Here, aaioo) = 
a'^{(jj)^ia'l^{uj) is acceptor's complex dipole polarizability, 
Qpiuo) and Rpiuj) are GP wavenumber and characteristic 
travel length, respectively; fdij^) is donor's normalized 
emission spectral function, and k is the effective dielectric 
constant of the environment (/^ = 1 for vacuum and k = 
2.5 for Si02 substrate). 

Figure [2] shows our numerical and analytical results for 
ETE between PbSe QDs a near graphene sheet doped to 
the Fermi level of e/ = 0.6 eV for several values of elec- 
tron scattering rate 7. Numerical results are obtained us- 
ing the full graphene density correlation function, while 
analytical results, given by Eqs. ([T]) and ([2|, are obtained 
within the plasmon pole approximation (see Section Il|) ; 
they are in excellent agreement for distances exceeding 
GP wavelength \p ^ 30 nm. The large-distance behavior 
of ETE depends strongly on the sample quality charac- 
terized by 7, which, in turn, determines GP travel length 
Rp. As is seen, for low-7 samples, GP-assisted ETE ex- 
ceeds Forster 's ETE in vacuum with ri? = 8 nm (shown 



10" 



10"' 



10 



10 



o 7 = 0.5 ps ^ 
□ 7 = 10 ps"^ 
A 7 = 40 ps~^ 



0.001 




_l I I I I I I I 



J I I I I I I I I 



0.01 ^ 0.1 

Rad (/^m) 



FIG. 2. ETE versus in-plane distance between donor and 
acceptor near suspended graphene in vacuum {k = 1). The 
Fermi level is set to e/ = 0.6 eV. Analytical and numerical 
results are shown by symbols and solid lines, respectively. 
ETE with and without graphene is shown in color and black, 
respectively. Further details are provided in the text. 



by the solid black line) by a factor of ~10^ for Ra 
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nm and by factor of ~10^ - 10^ for Rad = 500 nm. Such 
extreme enhancements are due to the slow decay of two- 
dimensional GP wave amplitude. For high-7 samples, 
ET is limited due to reduced GP travel length Rp. 

In the rest of the paper, we derive Eqs. ([T]) and ([2|, and 
study numerically and analytically ET over graphene for 
various system parameters. 



II. THEORY 

We consider donor and acceptor fluorophores (dye 
molecules or QDs) as point dipoles situated at = 
(Ri,Zi), {i = a,(i) with transition dipole moments 

= /iiiii, where is the dipole orientation, sepa- 

1/2 

rated from each other by Tad = [^ad + i^a ~ ^d)'^] 
(see inset in Figure 2|. If 2:^ is not too small so that 
fluorophores' internal transitions are not signiflcantly af- 
fected by graphene, ETE can be found within the semi- 
classical approach. The power transferred from the 
donor, initially excited at frequency cj, to the acceptor is 
given by 
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Pad{^) = 2 ^a(^) |na ' E(a;; r«) 



(3) 



where E(cc;; r^) is the electric fleld at the acceptor's posi- 
tion. This fleld is related to the donor's dipole moment 
via 



E(cj; r) = — 5— G(cj; r, r^) • /j.d, 



(4) 



where G(co'; r, r') = Go{uj; r, r') + Gg{oo; r, r') is the elec- 
tric fleld Green dyadic comprised of direct and graphene- 
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assisted contributions, respectively. For brevity, we intro- 
duce a matrix Sij{uj) = (47ra;^/c^) • G{uj;Ti,Tj) • n^, 
with similar decomposition Sij = Sfj{uj) + Sfjiuo). In 
terms of Sij^ the transferred power, Eq. (|3|, takes a sim- 
ple form Pad = {uj/2)iila'^ \Sad\^ - 

ETE is obtained by normalizing Pad{^) with respect to 
the donor's full power loss, Pd{^), followed by integration 
over the donor's emission band: Ead = f dujfdPad/Pd- In 
the lowest order, Pd{oj) has the form 



P4 = p', 
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where P^ stands for the donor's power loss due to ra- 
diative and intrinsic (non-radiative) processes and P^ is 
the power dissipated in graphene. In vacuum, the for- 
mer is given by P^ = {uj/2Qd)fillmS^^ = fijuj^ /Sc^Qd, 
where the donor's quantum yield, Qd^ accounts for in- 
trinsic losses, while in the presence of dielectric inter- 
face (formed by an under-graphene substrate, e.g., Si02), 
it is more involvec^ and has been evaluated by us nu- 
merically. The power dissipated in graphene is given by 
P| = (a;/2)/i^Im6'^^. For Zi > 1 nm considered here, 
higher order terms describing feedback from acceptor to 
graphene and from graphene to donoi^ are small and, 
therefore, neglected. The ETE then takes the form 
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where averaging over dipoles' orientations is implied (see 
Appendix [B] for the detailed discussion of the ETE de- 
pendence on specific dipoles' orientations). 

We now proceed with evaluation of 5'^^ = S^^^uo) + 
Sf-{uj). The direct (Forster) ETE is determined by the 

Coulomb term in free space, S^^ = Qad/'^ad^^^ where Qad is 
the orientational factor with average (g^^) =2/3, while 
the donor's radiative losses are described by Im^^^ = 
^{uj/c)^. The graphene contribution to Sij can be found 
as follows. In the longwave limit where retardation effects 
can be neglected, Sad{^) reduces to 
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where U = v-\- vHv is the Coulomb potential screened by 
the graphene sheet, v{r) = v{Il^ z) is the bare Coulomb 
potential and n(R, cj) is the density correlation function 
of graphene. After the in-plane Fourier transform using 

as 



l^e ^1^1, the graphene contribution is obtained 
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where ^i(q) = • q + • is the orientational fac- 
tor, q and z^ being, respectively, the radial unit vector 
in the graphene's plane and the normal from graphene's 
plane to fiuorophore i. Eqs (|6| and ^ are used by us to 
numerically evaluate ETE and obtain all the numerical 
results in this work. Specifically, all the Green dyadics 



in the matrix representation {Sadi S^d ^dd) ^^^^ 
evaluated (for each uj) via numerical integration over the 
wavenumber q adopting the density correlation function 
of the homogeneous graphene, 11(^^,0;), in the random 
phase approximation (see Appendix [A| for details). Then, 
the integration over uj in Eq. (6| is performed numerically. 

Analytical expression for the long-distance behavior of 
5*^^ can be derived using the plasmon pole approximation 
for n(g,cj) as (see Appendix [A|) 



n(^,c^)^n^^(g,a;) 
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where Ap, Qp and Rp^ being respectively the CP ampli- 
tude, wavenumber and travel length, are obtained by lo- 
cating the resonance of 11(^,0;) at g = Qpi^)- At low 
energies (i.e., huj <C e/), this can be done analytically 
yielding 
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qp = kh UJ /2e/e 



Rp = CfC^ /kh^ujj. 



(10a) 



(10b) 



(10c) 



It turns out, that even at ficj ~ e/ (i.e., regime consid- 
ered in this work) the low-cj analytical expressions for 
Ap and qp are still applicable. In contrast, Rp has to be 
found very accurately (i.e., numerically) since Eq. (|2| IS 
exponentially sensitive to its value. Our additional nu- 
merical tests (not shown) have demonstrated that the 
analytical expression for Rp becomes quite accurate al- 
ready at hw/ef < 0.3 - 0.5. However, at fkj ~ ej one 
can expect the magnitude of the inaccuracy of Rp^ as ob- 
tained from Eq. ( |10c[ ), to be of the order of the value of 
Rp. Specifically, for the parameters adopted in this pa- 
per, the analytically found Rp is approximately twice as 
high as its numerical counterpart. In what follows, the 
analytical expressions for Ap, qp^ and numerically found 
Rp are adopted. 

Substitution of Eq. (|9| into Eq. (|8| yields 



dqq 



^-q{\Za\ + \Zd\) 

Qp - i/2Rp 



d(l)fa{^)fd{^)e 



iq R 



(11) 

where (j) is the azimuthal angle in q-plane. For qR^ 1^ 
only small fiuctuations of (j) around q • R = con- 
tribute to ^-integral /, and in these regions /i(q) can be 
replaced by /i(±R), yielding 
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Upon substitution of Eq. (12) into Eq. (11), S^^ splits 



into two parts corresponding to the outgoing and incom- 
ing waves, e^*^^. For qR ^ 1, the dominant contribution 
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comes from the pole at g = + i/2Rp into the outgoing 
part. Finally, after averaging of \S^^J^ over dipoles' ori- 
entations using relations (fifj) = §Sij and (fifj) = 0, 
one obtains 
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Comparison of Eq. (13) and direct contribution l^^^ 



i^ad reveals that the GP-assisted ET channel is dom- 
inant for Rad ^ Ap. Specifically, numerical calcula- 
tions point to a crossover to the GP-assisted regime at 
Rad -- 30 - 40 nm for ^ 30 nm (see Figure ||). GP- 
assisted ET is ineffective for large fiuorophores' distances 
to the graphene plane {ziQp ^ 1) or for their in-plane 
separation significantly exceeding plasmon travel length 

{Rad/Rp > 1). 

Turning to dissipated power in graphene, Pj', the di- 
agonal element S^^ can be obtained from Eqs. ([8| and 
(|9| via substitution a ^ which yields Im6^^ = 

^^^-'^qp\zd\S5 j£ acceptor is absent, P| determines 
the normalized energy loss rate of the donor, kd = 
(P^ + Pj'j/P^, shown in Figure [IJc) for several values 
of hu. At very large Zd, the donor does not "feel" the 
presence of graphene, so its losses are dominated by ra- 
diative and non-radiative ones. At smaller distances, the 
exponential decay of kd with the donor-graphene distance 
indicates the predominant donor's energy transfer to GP. 
At even smaller fiuorophore- graphene distances, the non- 
exponential dependence of kd on Zd is due to the onset of 
excitations of electron-hole pairs in graphene [see inset in 
Figure [ijc)]. These three regimes of a single fiuorophore 
interaction with graphe ne have recently been studied in 
detail elsewhere .ESSEH 

The above considerations lead to a conclusion that in 
the wide range of intermediate donor-graphene distances 
(i) donor's energy losses are dominated by GP excitation, 
and (ii) kd ^ 1 and so P^ dominates over intrinsic and 
radiative losses. Furthermore, at distances between fiu- 
orophores exceeding Forster radius, dominates over 
Pad in Eq. ([5| as well. Thus, in a wide parameter range, 
both the numerator and denominator of the integrand in 
Eq. ([gI are dominated by GP-assisted channels, yielding 
Eqs.^ and ^ for ETE. 



III. RESULTS AND DISCUSSION 



ET calculations below were performed for a pair (donor 
and acceptor) of PbSe QDs with emission and absorption 
bands centered at 0.55 eV and 0.6 eV, respectiveljS^ED. 
The fiuorescence quantum yield for such QDs varies sig- 
nificantly in literature, so the "average" value of 10~^ 
is adopted here. Lorentzian lineshape for both bands 
is assume d with full width at half maximum (FWHM) 
of 0.1 eV j^^*^^ * ^ ^ and the acceptor absorption crossec- 
tion is ch osen ctq = (Attuj/Sc) a^^ = 2 A^ at its spectral 
maximum .E^'^^EH Both optical bands lie within the GP 



band with dispersion uu ex ^ in doped graphene with 
electron scattering rate chosen as 7 = 10 pg-i 132J361 Por 
the Fermi level at e/ = 0.6 eV adopted here,EHlGP is well 
defined up to q ~ 0.6 nm~^ corresponding to hu ~ 0.8 
eV, while for larger q GP is dampened by interband 
single-particle transitions - Landau damping [see Fig- 
ure fljb)]. A donor with emission band centered at 0.55 
eV [dashed horizontal line in Figure ijb)] predominantly 
excites CPs with qp ^ 0.2 nm~^ (red dagger), while exci- 
tation of electron-hole pairs requires higher wavenumbers 
g > 0.9 nm~^ (blue half-oval) and is, therefore, efficient 
only for Zd < I nm [see inset in Figure [ijc)]. Below 
we choose the values Zd = Za = 3 nm lying in the GP- 
dominated exponential domain with kd ^ ^ x 10^ [see 
Figure [l{c)]. 

In Figure [2j the results of our numerical and analyti- 
cal calculations of ETE for suspended graphene are com- 
pared to Forster's ETE for a similar system in vacuum. 
Forster's ETE shows characteristic behavior described by 
standard expression E^^ = (l + r^d/^r) ' where calcu- 
lated Forster radius r^? ^ 8 A is consistent with experi- 
mental results for a similar system.'^ Numerical results for 
ETE in the presence of graphene (solid lines) are shown 
for several values of electron scattering rate 7. For small 
Rad^ the energy transfer from donor to acceptor, deter- 
mined by the integrand's numerator in Eq. (6|, is dom- 
inated by the direct Forster mechanism. However, the 
donor energy losses^ defined by the integrand's denomi- 
nator, are greatly increased, as compared to the vacuum 
case, due to the presence of graphene. Under these condi- 
tions, the ETE dependence on the distance between QDs 

becomes Forster-like again with Ead = (l + ^ad/^^) 
However, the effective transfer radius, ~ 2 nm, is now 
significantly smaller than rp in the vacuum case due to a 
much larger, compared to radiative and intrinsic losses, 
donor energy dissipation to graphene. 

For large Rad^ ETE exhibits significant dispersion for 
different values of 7 caused by reduction of plasmon travel 
length Rp with increasing 7 [see Eq. (I9|] and, hence, the 
exponential suppression of ETE for Rad ^ Rp- The dis- 
tance dependence of ETE for Rad ^ 10 nm is in excellent 
agreement with our analytical results, Eqs. ([T]) and ([2|, 
shown by symbols in Figure [2j 

In Figure [3] we show how ETE evolves with the dop- 
ing level of graphene. As the Fermi level is reduced from 
0.6 eV to 0.2 eV with the decrement of 0.1 eV, ETE 
first decreases slowly and then sharply drops at e j below 
0.5 eV, i.e., when GP Landau damping onset (~1.3e/) 
moves below donor's emission band. For Cf = 0.2 eV, 
ETE shows Forster-like behavior (xR~^ but with the re- 
duced effect ive radiu s of ^2.5 nm due to ET quenching 
by graphene.l^^ESIlSI 

In Figure [4j we plot calculated ETE vs. fiuorophores' 
separation from the graphene sheet (za^Zd) for different 
values of in-plane distance Rad- Here, we distinguish be- 
tween two scenarios: (a) the acceptor's position is fixed, 
whereas the donor is moved away from the graphene 
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FIG. 3. Dependence of ETE between two PbSe QDs situated 
near graphene on the doping level of graphene. The adopted 
parameters are 7 = 10 ps~^ and Zd = Za = S nm. The 
schematic of the system is shown in the inset. 
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FIG. 4. Dependence of ETE on the distance between 
donor/acceptor and graphene. Legend encodes Rad- (a) 
Donor-graphene distance, Zd, is varied, while the acceptor- 
graphene distance kept constant (za = 3 nm). (b) Both dis- 
tances are varied simultaneously {za = Zd). Doped graphene 
(e/ = 0.6 eV) lays on top of the Si02 substrate {k = 2.5). 
The adopted electron scattering rate is 7 = 10 ps~^. The 
gradient-colored bar between the panels marks the transition 
from the GP-mediated ET at small Zd to the standard FRET 
mechanism at larger Zd- 



plane, so that donor-acceptor distance Tad increases as 
well; and (b) the acceptor follows the donor so that both 
fluorophores are moved in sine away from the graphene 
plane, i.e., rad stays constant for each Rad value. In both 
scenarios, at Zd ^ 10 - 20 nm both ET and losses are 
dominated by the GP-assisted channel, and, therefore. 



Eqs. ([T]) and ([2| are expected to provide an accurate de- 
scription of ETE behavior. Indeed, ETE plateaus [panel 
(a)] and the exponential decay of ETE [panel (b)] at low 
Zd both originate from the same exponent in Eq. ([2|, 
which is independent of Zd and linear with respect z^, 
respectively. 

At large Zd, the GP amplitude is exponentially damp- 
ened [Eq. ([8|], i.e., graphene becomes effectively absent 
from the ET picture, so that ETE dependence on Zd fol- 
lows the standard FRET-like behavior. Specifically, 
Ead R~^ behavior for Za = Zd results in plateaus with 
i^ad-dependent levels at large Zd [panel (b)]. In panel (a), 
this dependence reduces to z^^ at very large Zd (i.e., 
Zd ^ Rad^ ^a)' If the GP-assisted channel is already neg- 
ligible but Zd is still much smaller than Rad - this regime 
can be realized at Rad ^ 30 nm - then Vad ^ Rad 
and Ead levels off with respect to Zd at, e.g., = 30 - 
100 nm for Rad = 100 and 200 nm [panel (a)]. Note that 
the magnitudes of large- 2:^^ plateaus in panels (a) and (b) 
match for each Rad value. 

Finally, Figure [4]^a) shows that the transition from 
GP-dominated to FRET-dominated ET results in the 
ETE increase for Rad = 10-20 nm and its decrease 
for larger in-plane distances. Bearing in mind the ef- 
fective "absence" of graphene at large Zd, this behav- 
ior can be traced back to that in Figure [2j where ET 
without graphene (i.e., FRET) is more efficient than 
the GP-mediated ET in the presence of graphene at 
Rad ^ 30 nm, and less efficient for larger in-plane dis- 
tances. 



IV. CONCLUSION 

In this paper, we have shown that a single- atom layer 
of doped graphene can be used for highly efficient long- 
range energy transfer at the nanoscale. The transfer 
is mediated by longitudinal plasmons in graphene and 
hence it is very sensitive to the sample mobility and dop- 
ing level which determine plasmon lifetime and travel 
length. We have demonstrated that in clean samples with 
high doping levels (e.g., = 0.6 eV), the energy transfer 
efficiency can exceed that of FRET by up to ^10^ at hun- 
dreds nm distances. For a given donor-acceptor pair, the 
transfer efficiency can be optimized by tuning parameters 
of the system, e.g., fluorophore-graphene distances. 
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Appendix A: Density correlation function 

The bare density correlation function, or retarded po- 
larization operator, is calculated within the Dirac elec- 
trons approximation a J^*^-*^'^^' 



Scf G(-A_)6>[-Re{A_}-l] 



2^2 



hvjq 



[G(A_) + iir] [Re {A_} + 1] - G(A+) 



(Al) 



where G{z) = z\J — 1 — In + — l) and A± = 
{uj/vf ^2ef /hvf) /q. The square roots are chosen to 
yield positive real parts and the imaginary part of the 
logarithm is taken in (— 7r,7r] range. Fermi velocity and 
Fermi level (the latter determines the extent of graphene 
doping) are denoted hy Vf and ej, respectively. Within 
the Dirac electrons approximation, the density correla- 
tion function is insensitive to the sign of the Fermi level, 
so in all the expressions here and in the main text e/ has 
to be understood as |ej|. 

The two important limiting forms of the density corre- 
lation function are (i) the long wavelength limit ^ 0, 
hjo <C 2e/), and (ii) the static limit {uo 0, q < 2k j). 
The long wavelength limit is given by 



Uo{q ^ 0,cj) = 



(A2) 



The static limit of the bare density correlation function 
is obtained as 



Uo{q,uJ 0) 



2e 



/ 



(A3) 



The naive substitution uj ^ uj-\- ^7/2 to account for in- 
graphene scattering losses in Eq. ( |A1[ ) (7 is the electron 
scattering rate) is inaccurate in a general case (especially 
if 7 is not small), since it does not preserve the particle 
conservation requirement. To correct for this, the m ore 
accurate Mermin procedure is adopted, yieldin^^^^ 



U^{q,uj) 



(1 + ij/uj)Uo{q,uj ^ ij) 
1 + {ij/u)Uo{q, u + ij)/Uo{q, 0) * 



(A4) 



The full (or "dressed") density correlation function, 
which accounts for screening in graphene, is obtained 
within the random phase approximation as 



U{q,w) = 



1 - e'^v{q)n^{q,uj)' 



(A5) 



where v{q) =2ti /kq is the two-dimensional Fourier trans- 
form of the Coulomb potential within the graphene's 
plane, v{R) = l/kR. The effective dielectric constant 
of the environment is given hy R = {hci -\- f<i2)/2 for a 



graphene sheet sandwiched between two homogeneous 
dielectrics with dielectric constants hzi and /^2-^^ Thus 
= 1 for a suspended graphene sheet in vacuum. For 
graphene, laid on top of a Si02 substrate (a^i = 1, 
^2 = ^Si02 = 4), one obtains R = 2.5. 

The plasmon dispersion relation, qp = qp{uo), is found 



by requiring the real part of the denominator of Eq. ( A5 ) 
to vanish. The Taylor expansion of the denominator 
around this point (up to leading terms in both real and 
imaginary parts) leads to the possibility of approximating 
the full density correlation function within the so called 
plasmon pole approximation as 



A„ 



q + i/2Rp 



(A6) 



where Ap = Ii^{qp^uj) /A is the plasmon amplitude, and 
Rp = A/2B is the plasmon travel length. The coefficients 



of the Taylor expansion of the denominator of Eq. ( A5 ) 



are 



dq 



{v{q)Re [n^(g,c^)]) 



(i=(ip 



B 



^v{qp)lm \U.^{qp,u)\ . 



(A7) 



In the low-cj limit (i.e., hw <C e/), the density correlation 
function in the plasmon pole approximation can be ob- 
taine d pu rely analytically by (i) substituting Eqs. (A2) 
and (A3) into Eq. (A4), and (ii) using the so obtained 
n^(g, cj) to evaluate the Taylor expansion coefficients A 
and B. The result of these manipulations - the explicit 
analytical formulas for Ap, qp and Rp - is provided in the 
main text after Eq. ([9|. 

It turns out that for the specific case considered here, 
i.e., the plasmon pole approximation in the long wave- 
length limit, the same analytical expression for Ii^^{q^uj) 
could have been obtained in the limit of small 7 by us- 
ing the substitution uo u -\- ij/2 instead of the more 
general Mermin's procedure. It has to be emphasized, 
however, that such an agreement is not general and hard 
to foresee. Therefore, the more accurate Mermin's proce- 
dure has to be favored over more approximate methods 
of introducing the finite scattering rate into the density 
correlation function.™ 



Appendix B: Dipoles' orientation dependence of 
ETE 



Figure |5] shows the distance dependence of ETE for flu- 
orophores near graphene at different dipole orientations. 
Graphene is laid on top of the Si02 substrate {k = 2.5) 
and the Fermi level is set to ej = 0.6 eV. The electron 
scattering rate is assumed 7 = 10 ps~^. 

Dash-dotted and dash-double-dotted lines represent 
the both donor and acceptor transition dipoles fixed in 
the x-direction [schematically shown in Figure [TJa) in 
the main text] and in the z-direction (dipoles perpendic- 
ular to the graphene's plane), respectively. ETE in the 
projection-averaged case is shown by solid lines. 
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FIG. 5. ETE without (black) and with graphene (color lines) 
laid on top of the Si02 substrate {k ~ 2.5). The adopted pa- 
rameters are e/ — 0.6 eV, 7 = 10 ps~^ and Za — Zd =3 nm. 
Dash-dotted, dash-double-dotted, and solid lines correspond 
to cases where both the donor and acceptor transition dipole 
vectors are elongated in x- and 2;-directions, or averaged over 
all the directions, respectively. Dashed lines represent sum- 
ming over the acceptor's dipole projections instead of averag- 
ing. 



Dashed lines represent the case where the summation 
is performed over the acceptor's dipole projections in- 
stead of averaging (solid lines). This has to be done if 
the acceptor's dipole polarizability is isotropic, which is 
frequently the case for spherical semiconductor QDs as 
fluorophores. For example, this is true for PbSe QDs, 
where the dipole polarizability is isotropic due to the 
presence of four degenerate L-valleys corresponding to 



the four equivalent (111) directions in the face-centered 
cubic lattice of lead chalcogenides.l^ Within the analyt- 
ical plasmon pole approximation, to substitute the aver- 
aging over acceptor's projections with summation it suf- 
fices to add an extra factor of 3 into Eq. ( 13 ) in the main 
text. 

For comparison, the dependence of the Forster ETE 
(i.e., in the absence of graphene) on dipole orientations 
is shown by black lines. In this case, the standard Forster 
ETE with Ead = is recovered. The slightly 

smaller Forster radius, rp ^ 7.5 nm for the solid black 
line, than the one obtained for the suspended graphene 
{tf ^ S nm) (Figure [2] in the main text), is due to the 
Si02-induced dielectric screening = 2.5). 

Specific dipole orientations can lead to strong ETE 
variations in the crossover region between Forster- 
dominated and GP-dominated regimes. This behavior 
is due to the interference between the Forster and GP 
contributions to ET in the region where the magnitudes 
of these two contributions are comparable. In particular, 
the negative and positive interferences are seen for dipole 
projections fixed in z and x-directions, respectively, at 
Rad ^ 10 nm in Figure |5j 

At large donor-acceptor distances, a specific dipole ori- 
entation has no significant effect on ETE, except for the 
overall numerical factor of the order of ~1. For example, 
the donor with the transition dipole fixed in 2:-direction is 
twice as efficient in exciting GP than that with the dipole 
in x-direction.^^ However, the z-dipole excites plasmons 
isotropically within the graphene plane, whereas the GP 
emission of x-dipole has a characteristic dipolar pattern 
[see Figure [TJ a) in the main text], concentrated in the 
direction of the acceptor (and also in the opposite di- 
rection). This results in the same power transfer in both 
cases, but with lower power losses in the case of the dipole 
fixed in the x-direction, which ultimately yields twice as 
high ETE for x-dipoles than for z-dipoles. 
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